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Abstract 

Activation  and  concentration  polarization  effects  in  anode- supported  solid  oxide  fuel  cells  (SOFC)  were  examined.  The 
anode  and  the  cathode  consisted  respectively  of  porous,  composite,  contiguous  mixtures  of  Ni  + yttria-stabilized  zirconia 
(YSZ)  and  Sr-doped  LaMn03  (LSM)  + YSZ.  The  composite  electrode  provides  parallel  paths  for  oxygen  ions  (through 
YSZ),  electrons  (through  the  electronic  conductor;  Ni  for  the  anode  and  LSM  for  the  cathode),  and  gaseous  species  (through 
the  pores)  and  thereby  substantially  decreases  the  activation  polarization.  The  composite  electrode  effectively  spreads  the 
charge  transfer  reaction  from  the  electrolyte /electrode  interface  into  the  electrode.  At  low  current  densities  where  the 
activation  polarization  can  be  approximated  as  being  ohmic,  an  effective  charge  transfer  resistance,  is  defined  in  terms 
of  various  parameters,  including  the  intrinsic  charge  transfer  resistance,  Rcl,  which  is  a characteristic  of  the  electrocatalyst/ 
electrolyte  pair  (e.g.  LSM/YSZ),  and  the  electrode  thickness.  It  is  shown  that  the  /?'”  attains  an  asymptotic  value  at  large 
electrode  thicknesses.  The  limiting  value  of  /?'"  can  be  either  lower  or  higher  than  Rcl  depending  upon  the  magnitudes  of  the 
ionic  conductivity,  <rn  of  the  composite  electrode,  the  intrinsic  charge  transfer  resistance,  Rct,  and  the  grain  size  of  the 
electrode.  For  an  Rcl  of  1.2  ilcm2,  cr  of  0.02  S/cm  and  an  electrode  grain  size  of  2 |xm,  the  limiting  value  of  is  0.14 
ficm2  indicating  almost  an  order  of  magnitude  decrease  in  activation  polarization.  The  experimental  measurements  on  the 
cell  resistance  of  anode-supported  cells  as  a function  of  the  cathode  thickness  are  in  accord  with  the  theoretical  model.  The 
concentration  polarization  is  analyzed  by  taking  into  account  gas  transport  through  porous  electrodes.  It  is  shown  that  the 
voltage,  V,  vs.  current  density,  i,  traces  should  be  nonlinear  and  in  anode-supported  cells,  the  initial  concave  up  curvature 
(d2V/d i2  SO)  has  its  origin  in  both  activation  and  concentration  polarizations.  The  experimental  results  are  consistent  with 
the  theoretical  model.  © 2000  Elsevier  Science  B.V.  All  rights  reserved. 
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1.  Introduction 

Solid  state  devices  such  as  solid  oxide  fuel  cells 
(SOFC)  consist  of  a cathode,  an  electrolyte,  and  an 
anode.  Two  basic  designs  have  been  explored  in  the 
development  of  the  SOFC;  the  electrolyte-supported 
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and  the  electrode-supported  [1],  In  the  former,  the 
electrolyte  is  the  thickest  component  and  is  effective- 
ly the  support  structure.  In  electrolyte-supported 
cells,  the  thickness  of  the  electrolyte,  typically  YSZ. 
is  150  |xm  with  thin  electrodes  screen-printed  on  it 
[2,3],  In  the  latter,  one  of  the  two  electrodes,  either 
the  cathode  or  the  anode,  is  the  thickest  component 
and  the  support  structure. 

In  electrolyte-supported  cells  the  ohmic  contribu- 
tion is  large  due  to  high  electrolyte  resistivity.  For 
this  reason,  such  cells  are  being  developed  for 
operation  at  ~1000°C  where  the  electrolyte  resistiv- 
ity is  low,  typically  ~20  Hem.  In  cathode-supported 
cells,  a 30-40  |xm  layer  of  YSZ  is  deposited  on  a 
porous  LSM  cathode  of  ~2  mm  thickness  [4]. 
Similarly,  in  anode-supported  cells,  a YSZ  layer  of 
10-20  |xm  thickness  is  deposited  on  a relatively 
thick  Ni  + YSZ  anode  [1,5], 

The  overall  performance  of  such  devices  is  dic- 
tated by  various  polarizations;  namely  ohmic,  activa- 
tion, and  concentration.  The  main  contribution  to  the 
ohmic  polarization  is  due  to  the  electrolyte.  A high 
ionic  conductivity  and  a small  electrolyte  thickness 
are  the  desired  characteristics  of  the  solid  electrolyte 
to  minimize  the  ohmic  contribution.  Although  vari- 
ous solid  electrolytes  with  high  ionic  conductivities 
at  moderate  temperatures  have  been  explored,  yttria- 
stabilized  zirconia  (YSZ)  is  by  far  the  most  widely 
used  solid  electrolyte  due  to  its  excellent  stability  in 
both  reducing  and  oxidizing  environments,  even 
though  its  conductivity  is  lower  (resistivity  is  higher) 
than  other  materials  such  as  ceria  and  Sr-  and  Mg- 
doped  LaGaO,  (LSGM).  In  electrolyte-supported 
cells,  the  typical  YSZ  membrane  thickness  is  150 
|xm.  At  800°C,  the  resistivity  of  YSZ  is  about  50 
flcm  which  translates  into  an  area  specific  resistance 
(electrolyte  contribution)  of  about  0.75  Ocm2.  This 
value  is  very  high  with  the  result  that  the  power 
density  is  rather  low  at  temperatures  below  about 
950°C.  However,  in  electrode-supported  cells,  the 
YSZ  electrolyte  thickness  need  only  be  about  10  pm. 
With  a 10  pm  thick  YSZ  film  as  the  electrolyte,  the 
area  specific  resistance  of  the  electrolyte  is  only 
about  0.05  Hem2  at  800°C  thus  making  an  efficient 
operation  of  an  SOFC  at  such  a low  temperature  in 
principle  possible.  However,  it  is  usually  observed 
that  despite  a low  ohmic  contribution,  the  area 


specific  resistance  of  the  cell  as  a whole  may  be 
several  times  larger.  The  reason  is  that  activation  and 
concentration  polarizations  can  often  outweigh  the 
ohmic  contribution.  Thus,  the  minimization  of  the 
overall  area  specific  resistance  of  cells  requires  an 
optimization  of  parameters  which  govern  the  activa- 
tion and  concentration  polarization  effects. 

Activation  polarization  is  related  to  charge  transfer 
processes  and  thus  depends  upon  the  nature  of 
electrode-electrolyte  interfaces.  In  composite  elec- 
trodes comprising  a mixture  of  an  ionic  conductor 
and  an  electronic  conductor,  or  in  a single  phase 
mixed  ionic  electronically  conducting  (MIEC)  elec- 
trode, the  process  of  charge  transfer  is  expected  to 
occur  over  some  distance  from  the  electrolyte/elec- 
trode interface  into  the  electrode.  That  is,  effectively 
the  interface  is  diffuse  insofar  as  the  region  over 
which  charge  transfer  occurs.  Concentration  polari- 
zation is  related  to  the  transport  of  gaseous  species 
through  porous  electrodes,  and  thus  is  related  to  the 
microstructure  of  the  electrodes;  specifically  the 
volume  percent  porosity,  the  pore  size,  and  the 
tortuosity  factor.  It  is  thus  clear  that  a reduction  in 
the  overall  area  specific  resistance  of  cells  and  a 
concomitant  improvement  in  cell  performance  can  be 
achieved  only  through  a reduction  in  these  two 
polarization  effects.  Indeed,  recent  work  has  shown 
that  an  area  specific  resistance  less  than  0.15  fiem2 
can  be  realized  in  anode-supported  cells  at  800°C 
with  a maximum  power  density  in  excess  of  1.8 
fi/cm"  [6].  This  was  achieved  by  minimizing  both 
activation  and  concentration  polarizations.  Effective- 
ly, the  microstructures  of  the  electrodes  and  the 
electrolyte /electrode  interfaces  were  tailored  to  mini- 
mize these  polarization  effects.  In  this  paper,  these 
two  polarization  effects  are  examined  in  light  of 
anode-supported  SOFCs. 

2.  Activation  polarization 

The  role  of  composite  electrodes,  comprising  a 
mixture  of  a solid  electrolyte  and  an  electronic 
conductor,  or  a mixed  ionic-electronically  conduct- 
ing (MIEC)  material,  on  electrode  kinetics  has  been 
examined  by  a number  of  researchers  [7-11].  The 
rationale  for  the  use  of  a composite  or  an  MIEC 


A.V.  Virkar  et  at.  1 Solid  State  Ionics  131  (2000)  189-198 


191 


electrode  is  that  it  allows  a spreading  of  the  reaction 
zone  from  the  electrolyte /electrode  interface  into  the 
electrode.  That  is,  effectively  it  results  in  a sort  of  a 
diffuse  electrolyte /electrode  interface  wherein  the 
charge  transfer  reaction  occurs.  Even  though  it  has 
been  proposed  by  some  that  in  a mixed  conducting 
material  as  an  electrode  the  charge  transfer  reaction 
can  occur  over  the  entire  surface,  it  is  generally 
thought  to  occur  at  electrolyte-electrocatalyst-gas 
three  phase  boundaries  (TPB),  which  are  also  high 
energy  sites.  In  a single  phase  MIEC  electrode,  the 
charge  transfer  may  predominantly  occur  at  high 
energy  sites  such  as  grain  boundaries  and  defects  in 
the  MIEC  at  the  MIEC /gas  phase  interface.  Tanner 
et  al.  [11]  examined  charge  transfer  reactions  in  a 
composite  electrode.  For  the  case  of  SOFC,  for 
example,  the  analysis  is  applicable  to  Ni+YSZ 
anode  and  to  LSM  + YSZ  cathode.  It  is  known  that 
the  electrode  overpotential  at  a Ni  + YSZ  electrode  is 
usually  negligible  and  much  of  the  activation  polari- 
zation can  be  attributed  to  the  cathode.  In  this 
context,  the  following  discussion  is  particularly 
applicable  to  Sr-doped  LaMn03  (LSM)  + YSZ 
cathode.  The  electrocatalytic  properties  of  LSM  over 
YSZ  can  in  principle  be  described  in  terms  of  a 
phenomenological  model  such  as  the  Butler-Volmer 
equation  which  relates  the  current  density,  i,  to  the 
activation  overpotential,  by 

. . / ( azFrj^  ( - (1  - a)zFr)act\  \ 

1 = 1A^V~rt^J  - exPl rt ) J 

(1) 

where  i0  is  the  exchange  current  density,  a is  the 
transfer  coefficient,  z is  the  number  of  electrons 
participating  in  the  electrode  reaction,  F is  the 
Faraday  constant,  R is  the  gas  constant,  and  T is  the 
temperature.  In  the  low  current  density  regime, 


RT 

V*ct  = 


(2) 


where  Rct  given  by 
RT 


is  the  intrinsic  charge  transfer  resistance  in  Ocm’.  It 
is  a function  of  the  electrochemical  properties  of  the 
electrocatalyst /electrolyte  pair  (LSM /YSZ  here), 
and  also  a function  of  the  TPB  line  length.  Thus,  it  is 
a function  of  the  particle  size  of  LSM  and  the 
amount  of  LSM  per  unit  area  of  YSZ.  In  what 
follows,  Rct  will  be  treated  as  an  empirical  parame- 
ter, determined  experimentally  for  a given  electro- 
catalyst/electrolyte pair.  In  composite  electrodes 
comprising  a mixture  of  LSM  + YSZ  for  instance, 
the  reaction  zone,  that  is  the  region  over  which  the 
process  of  charge  transfer  occurs,  is  spread  out  from 
the  electrolyte /electrode  interface  into  the  electrode. 
In  such  a case,  it  is  expected  that  the  activation 
overpotential,  r;oct,  may  be  lower  than  that  given  by 
Eq.  (2).  Tanner  et  al.  [11]  defined  an  effective  charge 
transfer  resistance,  given  by 


^ac  = Klfi  (4) 

Using  this  as  a parameter  for  the  design  of  composite 
electrodes,  Tanner  et  al.  [11]  derived  an  expression 
for  Rec?  terms  of  Rct  and  the  ionic  conductivity  of  the 
electrode,  cr.,  which  describe  the  properties  of  the 
electrode /electrolyte  interface  and  the  electrolyte, 
respectively,  and  microstructural  parameters  of  the 
electrode,  namely  the  grain  size  and  porosity.  This 
equation  is  given  by  [11] 


Rt?  = ' 


BR„ 


(1  ~Vv)  exp 


(-!) 


1 + /?  exp 


+ 


1 + /?  exp 


Bl 

(-¥), 


a(i  — exp(— y))  +bV„ 
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where 


A =\](riB{l  — Vv)Rcl  and 


(t,Rci  ~ A 


(6) 


in  which  B is  the  grain  size  of  the  electrolyte  in  the 
composite  electrode,  Vv  denotes  the  fractional  po- 
rosity, and  h is  the  electrode  thickness.  It  is  to  be 
emphasized  that  the  ionic  conductivity  of  the  elec- 
trode, cr  in  general  may  be  different  from  that  of  the 
dense  electrolyte  membrane.  For  example,  the  elec- 
trolyte membrane  may  be  of  say  YSZ  but  the 
electrode  may  contain  doped  ceria.  The  analysis 
assumes  that  the  electronic  conductivity  of  the 
electrode  is  much  larger  than  the  ionic  conductivity, 
and  gas  transport  through  the  electrode  is  not  rate- 
limiting.  Eq.  (5)  shows  that  the  effective  charge 
transfer  resistance,  Rc"  approaches  /?cl  as  h -» 0,  and 
achieves  an  asymptotic  value  with  an  increasing 
electrode  thickness.  It  can  be  shown  that  as  h — > », 
the  /?*[f  is  given  by  [II] 

<7» 

as  long  as  a;  is  not  too  small.  As  <7;  — >0,  the 
asymptotic  value  of  /T',"  is  given  by 


eff.^ 

* K 


(8) 


If  ai  -»  0,  it  means  that  only  the  three  phase  bound- 
ary (TPB)  at  the  surface  of  the  electrocatalyst/dense 
electrolyte  membrane  (e.g.  LSM  and  the  dense  YSZ 
membrane)  contributes  to  the  charge  transfer  re- 
action; hence  the  presence  of  Vv  in  the  denominator 
of  Eq.  (8).  Eqs.  (6)  and  (7)  also  show  the  effect  of 
the  microstructural  dimension,  B,  which  is  essentially 
the  grain  size  of  YSZ  in  the  electrode.  As  an 
illustration,  let  us  choose  the  following  values:  Rct  = 
1.2  ficm2,  a;.  = 0.02  S/cm,  B = 2 pan,  and  =Vr  = 
0.35.  Then,  the  estimated  value  of  the  effective 
charge  transfer  resistance  ReJt'  is  0.14  Qcm2,  almost 
an  order  of  magnitude  decrease  in  the  charge  transfer 
resistance,  or  an  order  of  magnitude  increase  in  the 
effective  exchange  current  density,  in(eff).  This  shows 
the  profound  effect  of  composite,  mixed  conducting 
electrodes  on  electrode  kinetics. 

Eq.  (7)  shows  that  the  effect  of  composite  elec- 


trode need  not  always  be  an  enhancement  in  elec- 
trode kinetics.  Indeed,  Eq.  (7)  shows  that 


R\ 


if 


B 


<>;■(  i-K) 


However, 

R ..  ff 

— >/?>/?  if 

1/  — “cl  — “et  11 


B 

oid-V,) 


^c« 


That  is,  the  electrode  kinetics  could  actually  be 
suppressed  if  the  electrode  parameters  are  unfavor- 
able; e.g.  too  coarse  a microstructure  (large  B)  and/ 
or  too  low  an  ionic  conductivity,  a;,  of  the  elec- 
trolyte material  in  the  electrode  (or  too  low  an  ionic 
conductivity  of  a single  phase  MIEC  electrode).  As 
an  example,  a suppression  of  the  electrode  kinetics 
would  occur  if  an  SOFC  anode  is  made  of  Ni  + an 
inert  material,  such  as  a-alumina.  Since  a-alumina  is 
not  an  oxygen  ion  conductor,  us  cris  essentially 
zero.  In  such  a case,  the  electrodef-poianzation  should 
actually  increase  when  using  such  a composite 
electrode.  Fig.  1 shows  plots  of  /?'[f  vs.  log  h (Eq. 
(5))  for  various  values  of  crr  Note  that  the  Rc)'  can 
either  decrease  or  increase  depending  upon  the  value 
of  cr  for  given  values  of  Ra,  B and  V„.  The 
experimental  results  of  Kenjo  et  al.  |7J  on  composite 
electrodes  showed  that  in  some  cases,  the  polariza- 
tion resistance  decreased  with  an  increasing  electrode 
thickness,  while  in  other  cases  there  appeared  to  be 
no  effect  of  the  thickness.  These  observations  are  at 
variance  with  the  predictions  of  the  model  by  Kenjo 
et  al.  [7],  However,  their  results  are  entirely  con- 
sistent with  the  model  presented  by  Tanner  et  al. 
[11]. 

Eq.  (5)  shows  the  dependence  of  Rcct  on  h.  For 
values  of  h greater  than  the  critical  value,  /?', 
achieves  an  asymptotic  value  given  by  Eq.  (7).  An 
experimental  verification  of  the  dependence  of  the 
asymptotic  value  of  /?'['  on  electrode  thickness  was 
obtained  on  anode-supported  solid  oxide  fuel  cells 
with  a varying  cathode  thickness.  The  anode  (Ni  + 
YSZ)  thickness  was  typically  1 mm  with  an  elec- 
trolyte (YSZ)  thickness  of  ~20  p,m.  The  cathode 
thickness  was  varied  over  a range  between  ~0  to 
~85  pan.  The  cathode  was  prepared  by  first  forming 
a porous  layer  of  YSZ  on  the  dense  YSZ  electrolyte 
layer,  and  then  infiltrating  an  aqueous  solution  of 
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Fig.  1.  Theoretically  calculated  values  of  the  effective  charge 
transfer  resistance,  fl'f.Tor  an  intrinsic  charge  transfer  resistance, 
RcV  of  1.2  flcm2  for  several  values  of  the  ionic  conductivity,  a;, 
of  the  cornpc'M.teelectfdde.  At  high  values  of  <r(,  the  1?*"  decreases 
with  an  increasing  electrode  thickness,  h,  and  achieves  a limiting 
value.  However,  it  is  seen  that  for  low  values  of  a;.,  the  R‘c" 
actually  increases  with  the  electrode  thickness,  h.  That  is,  under 
such  conditions,  the  composite  electrode  may  actually  suppress 
electrode  kinetics  (enhance  electrode  polarization). 


La-nitrate,  Sr-nitrate,  and  Mn-nitrate.  After  salt 
infiltration,  the  cells  were  heated  to  1000°C  to 
decompose  the  salts  and  form  LSM.  The  voltage,  V, 
vs.  current  density,  i,  performance  curves  were 
measured  at  800°C  with  humidified  hydrogen  as  the 
fuel  and  air  as  the  oxidant.  Fig.  2 shows  the  voltage 
vs.  current  density  polarization  curves  for  a number 
of  cells  for  different  cathode  thicknesses.  It  is  clearly 
seen  that  the  performance  increases  with  an  increas- 
ing cathode  thickness  up  to  ~45  p.m.  For  the  cell 
with  a cathode  thickness  of  ~85  pm,  the  perform- 
ance is  somewhat  lower  than  for  the  cell  with  a ~45 
pm  cathode.  This  is  attributed  to  three  possible 
reasons:  (a)  The  method  of  deposition  of  the  porous 
YSZ  layer  of  the  cathode  was  not  optimized  which 
led  to  a highly  tortuous  path  for  the  transport  of 
gaseous  species.  This  presumably  led  to  a significant 
concentration  polarization,  especially  at  larger 
cathode  thicknesses,  (b)  In  this  nonoptimized  method 
for  the  preparation  of  the  cathode  layer,  some 
cracking  of  the  layer  with  cracks  parallel  to  the 
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Fig.  2.  The  experimentally  measured  voltage,  V,  vs.  current 
density,  i,  traces  for  anode-supported  cells  with  varying  cathode 
thicknesses.  Note  that,  in  general,  the  performance  increases  with 
an  increase  in  cathode  thickness. 


electrolyte /cathode  interface  is  possible.  This  may 
have  led  to  some  increase  in  the  resistance,  (c)  The 
method  of  introducing  the  salt  solutions  was  not 
optimized  which  may  have  led  to  a nonuniform 
distribution  of  the  subsequent  LSM  formed  after 
calcination.  It  is  nevertheless  clear  that  the  cell 
performance  indeed  increases  with  an  increasing 
cathode  thickness.  From  the  near  linear  regions  of 
the  voltage  vs.  current  density  traces,  the  cell  resist- 
ance was  measured.  From  this,  the  electrolyte  area 
specific  resistance,  //,.,  calculated  using  independent- 
ly measured  ai  and  the  electrolyte  thickness,  was 
subtracted.  Neglecting  concentration  polarization, 
this  gives  an  approximate  estimate  of  the  effective 
charge  transfer  resistance,  R^.  Fig.  3 shows  the 
estimated  (experimental)  vs.  the  cathode  thick- 
ness (on  a log  scale).  Also  shown  is  a plot  of  ReJt 
calculated  using  Eq.  (5)  with  Rcc=  1.2  flcm2,  o;  = 
0.02  S/cm,  B — 5 p,m,  and  1^  = 0. 3.  Note  that  the 
agreement  between  the  experimental  data  and  the 
model  is  good. 

The  preceding  discussion  assumes  the  activation 
polarization  to  be  ohmic;  i.e.  i7act  being  linearly 
proportional  to  the  current  density,  i,  which  is 
applicable  for  current  densities  less  than  the  ex- 
change current  density,  ia.  The  composite  electrode 
thus  increases  the  exchange  current  density,  which 
may  be  termed  the  effective  exchange  current  den- 
sity,  io(ef0,  given  by 
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Log|{Jh)  (h  in  pm) 

Fig.  3.  A plot  of  /?[ [f  vs.  cathode  thickness,  h (log  scale), 
calculated  using  ^,  = 1.2  flcm2,  er.  = 0.02  S/cm,  and  the  grain 
size  of  YSZ  in  the  cathode  =5  pirn.  Also  plotted  arc  the 
experimentally  estimated  (approximate,  because  the  V vs.  i traces 
are  nonlinear)  /?'"  from  the  data  in  Fig.  2.  The  data  point 
corresponding  to  zero  cathode  thickness  is  that  corresponding  to 
the  LSM  paste  directly  applied  on  the  dense  YSZ  layer,  without  a 
porous  YSZ  layer.  Since  the  data  point  corresponding  to  zero 
cathode  thickness  actually  represents  that  for  a small  but  a 
nonzero  thickness,  and  that  zero  thickness  can  not  be  shown  on  a 
log  scale,  this  data  point  is  shown  at  h = 1 pm. 


of  the  voltage  vs.  current  density  traces  of  cells. 
Thus,  it  is  tacitly  assumed  that  the  composite  elec- 
trode modifies  a and  b such  that  inlcrl  ) is  much  larger 
than  ia. 


3.  Concentration  polarization 

Fig.  4 shows  schematics  of  cathode-supported  and 
anode-supported  cells.  In  the  former,  the  cathode  is 
the  thickest  component  of  the  cell  and  in  the  latter 
the  anode  is  the  thickest  component.  Fig.  4 also 
shows  schematic  variations  of  the  partial  pressures  of 
the  various  gaseous  species.  The  analysis  of  con- 
centration polarization  thus  should  begin  with  the 
analysis  of  the  transport  of  gases  through  porous 
electrodes.  In  the  anode,  the  principal  gaseous 
species  are  H,  and  H;0,  and  in  the  cathode  they  are 
O,  and  N2.  The  basic  equations  describing  an 
isothermal  transport  of  gaseous  species  through 
porous  electrodes  for  a mixture  of  two  gases,  A and 
B,  are  given  by  [12,13] 

J,  = - D,VnA  + XASAJ  - (11) 


io(cfn  ZFRll(  (9) 

However,  at  large  current  densities,  either  the  Bu- 
tler-Volmer  or  the  Tafel  equation  will  be  applicable. 
Preliminary  calculations  have  shown  that  the  compo- 
site electrode  can  be  effective  in  decreasing  the  r7.iet 
in  a similar  manner,  although  an  explicit  analytical 
form  is  not  available.  In  the  Tafel  limit,  the  r/act  may 
be  given  by 

17., c = a + b In  i (10) 

where  the  parameters  a and  b are  influenced  by  the 
electrode  microstructure  and  the  thickness.  In  what 
follows,  it  will  be  assumed  that  the  activation 
polarization  is  adequately  described  by  the  Tafel 
equation  and  it  will  be  incorporated  into  the  analysis 


•^b  + XDSBJ  XnyH 


where 


(12) 


(a)  Anode -S u ppo rted  Cell  (b)  Calhodc-Supportcd  Cell 


Fig.  4.  Schematics  showing  the  variation  of  partial  pressures  of 
hydrogen.  pHA  and  water  vapor,  pHM.  through  the  anode,  and 
oxygen,  pOA  through  the  cathode  of:  (a)  an  anode-supported  cell, 
and  (b)  a cathode-supported  cell. 
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In  Eqs.  (11)  and  (12),  the  fluxes  of  A and  B are 
given  by  JA  and  JB,  respectively,  DAB(eff)  is  the 
effective  binary  diffusivity  (taking  into  account  the 
porosity  and  the  tortuosity  factor),  DAK(eff)  & Z>BK(eff) 
are  the  effective  Knudsen  diffusivities,  nA  & nB  are 
the  concentrations  of  A & B (#/cm3),  respectively, 
n = (nA  + nB),  XA  & Xa  are  the  mole  fractions,  Ba  is 
the  permeability  through  the  porous  electrode,  /x  is 
the  viscosity,  p is  the  total  pressure,  and  J is  the 
total  flux.  The  effective  diffusivities  are  given  by  the 
intrinsic  diffusivities  multiplied  by  the  porosity,  Vv, 
and  divided  by  the  tortuosity  factor,  r [12-14],  The 
parameters  <5A  & SB,  and  yA  & yB  are  as  defined 
above.  Eqs.  (11)  and  (12)  include  two  flux  contribu- 
tions; a diffusive  flux  and  a viscous  flux.  The 
diffusive  contribution  consists  of  two  terms;  a free 
molecule  or  Knudsen  flow  (the  terms  containing 
^AK(eff)  & ^BK(eff))’  an(i  a continuum  part  (the  terms 
containing  the  effective  binary  diffusivities).  The 
typical  total  pressure  on  the  anodic  and  the  cathodic 
sides  is  on  the  order  of  1 atm  or  greater,  and  the 
typical  pore  size  is  & 1 p,m.  Thus,  in  general  it  may 
be  assumed  that  0AK(eff)  & DBK(ef()  » DABieff). 
Thus,  it  will  be  assumed  that  DA  & DB  can  be 
replaced  by  DAB(eff)  [12],  That  is,  at  high  pressures, 
such  as  those  of  interest  here,  both  DA  & DB 
approach  DAB(eff),  the  effective  binary  diffusion 
coefficient,  and  SA  & <5B  — > 1,  and  yA  & yB  0.  For 
the  cathode,  the  DAB  is  the  binary  diffusion  coeffi- 
cient for  a mixture  of  oxygen,  02,  and  nitrogen,  N2. 
For  the  anode,  the  binary  diffusion  coefficient  is  that 
for  H2  and  H20.  At  room  temperature,  DQ^_N  ~ 
0.22  cm2/s  and  Oh2-h2o  0-91  cm2/s  [13].  If  the 
fuel  also  contains  CO  and  C02,  then  their  transport 


must  also  be  accounted  for.  The  dependence  of 
diffusion  coefficient  on  temperature  is  not  very 
strong  and  to  a first  approximation  it  may  be 
modified  using  the  kinetic  theory  of  gases.  This 
suggests  that  for  comparable  porosities  and  electrode 
thicknesses,  the  concentration  polarization  effects  in 
general  should  be  less  on  the  anode  side. 

If  the  electrode  porosity  and  microstructure  are 
functions  of  position,  and/or  /)B(eff)  is  composition- 
dependent,  the  variation  of  partial  pressures  of  the 
various  species  with  position  will  not  be  linear.  For 
this  reason,  the  schematics  given  in  Fig.  4 show 
nonlinear  variations  of  partial  pressures. 

In  a steady  state,  the  equality 

I)h2I  = I/h2oI  = 2l-/oJ  ~ 2 F (14) 


must  hold,  where  jH ^ jH  O and  )0  i are  respectively 
the  fluxes  of  hydrogen  through  the  porous  anode,  of 
water  vapor  through  the  porous  anode  and  of  oxygen 
through  the  porous  cathode,  and  NA  is  the  Avogadro 
number.  The  flux  equations  were  solved  for  the 
steady  state,  subject  to  the  above  flux  equality 
requirement  [6].  By  incorporating  the  activation 
polarization  (using  the  Tafel  equation)  and  the  con- 
centration polarization  (by  solving  the  flux  equations 
for  transport  through  porous  electrodes),  it  was 
shown  that  the  voltage  vs.  current  density  behavior 
can  be  adequately  described  by  [6] 


RT 


V(i)  = Ea  — iRi  — a — b In  i + -r=-  In 


RT 

+ 2p  In 


4^) 

P H20!as  / 


(15) 


where  E0  is  the  open  circuit  voltage  (V),  R:  is  the 
area  specific  resistance  (ohmic)  of  the  cell  (flcm2), 
p'H^  is  the  partial  pressure  of  hydrogen  in  the  fuel 
(outside  of  the  anode),  p°H^0  is  the  partial  pressure  of 
water  vapor  in  the  fuel  (outside  of  the  anode),  ias  is 
the  anode-limiting  current  density  (a  current  density 
at  which  the  partial  pressure  of  hydrogen  at  the 
anode /electrolyte  interface  is  nearly  zero),  Amp/ 
cm2,  and  ics  is  the  cathode-limiting  current  density  (a 
current  density  at  which  the  partial  pressure  of 
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Fig.  5.  The  voltage,  V,  vs.  current  density,  i,  trace  for  an  improved  anode-supported  cel!  tested  at  800°C  with  humidified  hydrogen  as  the 
fuel  and  air  as  the  oxidant.  The  anode  thickness  was  -1.04  mm  and  the  electrolyte  (YSZ)  thickness  was  ~10  pm.  The  symbols  are  the 
experimental  data  points  while  the  curve  is  the  best  fit  to  Eq.  (15).  The  parameters  corresponding  to  the  best  fit  are  given  in  Table  1. 


oxygen  at  the  cathode /electrolyte  interface  is  nearly 
zero),  Amp  /cm2.  The  ias  and  the  i'cs  can  be  given  in 
terms  of  the  respective  electrode  thicknesses,  the 
respective  effective  binary  diffusivities  (which  in- 
clude the  porosity  and  the  tortuosity  factors),  and  p°H , 
and  /?o2,  respectively,  where  p")t  is  the  partial 
pressure  of  oxygen  in  the  oxidant  (0.21  atm  for  air  as 
the  oxidant)  [6], 

The  data  given  in  Fig.  2 are  on  some  early  cells 
wherein  the  porous  YSZ  layer  for  the  prospective 
cathode  was  first  deposited  followed  by  the  intro- 
duction of  an  aqueous  salt  solution  of  the  cathode 
precursor.  The  microstructure  in  such  cathodes  was 
however  rather  coarse.  Eq.  (7)  shows  that  the  grain 
size  of  the  YSZ  in  the  cathode  should  be  as  small  as 
possible.  Thus,  in  later  cells,  a much  finer  cathode 
microstructure  was  developed.  Also,  in  these  cells, 
the  YSZ  electrolyte  was  typically  about  10  p.m  in 
thickness;  about  half  as  thick  as  in  cells  for  which 
the  performance  data  are  given  in  Fig.  2. 

Fig.  5 shows  a plot  of  V vs.  i for  an  improved 
anode-supported  cell  tested  at  800°C  with  humidified 
hydrogen  as  the  fuel  and  air  as  the  oxidant.  The 
electrolyte  thickness  was  about  10  pim  and  the  anode 
thickness  was  about  1 mm.  The  corresponding  power 
density  vs.  current  density  plot  is  shown  in  Fig.  6.  As 
seen  in  Fig.  6,  a power  density  as  high  as  1.9  fl/cm" 


was  measured.  Such  a high  power  density  was 
achieved  by  minimizing  both  activation  and  con- 
centration polarizations  through  an  optimization  of 
the  electrode  microstructure.  The  symbols  in  Fig.  5 
are  the  experimental  data  points  and  the  curve  is  the 
best  fit  to  Eq.  (15).  For  an  anode-supported  cell,  the 


Current  Density  (Amps/cm  2) 

Fig.  6.  The  power  density  vs.  current  density,  r,  trace  for  the  cell 
in  Fig.  5.  The  maximum  power  density  measured  is  ~1.9  fl/cirT. 
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anode  thickness,  Za,  is  much  greater  than  the  cathode 
thickness,  /c,  with  the  result  that  Zas  Zcs.  Thus,  the 
term  containing  Zcs  from  Eq.  (15)  can  be  neglected. 
Table  1 gives  the  parameters  corresponding  to  the 
best  fit  to  Eq.  (15),  as  well  as  some  of  the  ex- 
perimentally measured  parameters.  The  effective 
exchange  current  density,  io(ef0,  was  estimated  to  be 
-110  mA/cm2,  which  is  more  than  an  order  of 
magnitude  higher  than  that  corresponding  to  LSM  on 
YSZ  at  800°C.  The  electrolyte  area  specific  resist- 
ance, /?;,  from  the  best  fit  was  estimated  to  be 
-0.046  Hem2.  The  electrolyte  thickness  was  about 
10  (Jim.  With  the  resistivity  of  YSZ  as  -50  flcm  at 
800°C,  the  estimated  value  of  R:  is  about  -0.05 
flcm2,  in  good  agreement  with  the  value  obtained 
from  the  curve  fit.  Finally,  the  estimated  value  of  the 
effective  binary  diffusivity  on  the  anode  side,  Z)a(eff), 
is  about  0.31  cm2/s  This  value  is  reasonable  when 
the  porosity  and  the  tortuosity  factor  are  taken  into 
account. 

The  theoretical  analysis  predicts  that  the  initial 
concave-up  curvature  (d2V/di2  >0)  in  V vs.  i traces 
in  anode-supported  cells  is  due  to  two  reasons  [6]: 
(a)  An  activation  polarization  that  is  nonohmic  (e.g. 
Tafel).  (b)  A concentration  polarization  that  is 

Table  1 


Experimentally  measured  and  curve-fitted  parameters  for  the 
voltage  vs.  current  density  trace  shown  in  Fig.  5 


Parameter 

Parameter 

Measured 

From  the  curve 

(description) 

(symbol) 

fit  to  the  V vs.  i 

Anode  thickness 

ia 

1.04  mm 

Tafel  coefficient 

a 

0.115 

Tafel  coefficient 

b 

0.0552 

Effective  exchange 
Current  density 

^o(eff) 

110  mA/cm2 

Cell  area 
Specific 
Resistance 
(Ohmic) 
Short  circuit 

R. 

K 

6.35  A/cm2 

0.046  flcm2 

Current  density 
Anode  limited 

ias 

6.35  A/cm2 

Short  circuit 
Current  density 
Anode  effective 
Binary  diffusion 
Coefficient 
Maximum  power 
Density 

Aa(eff) 

1.9  W/cm2 

0.31  cm2/s 

associated  with  too  low  a partial  pressure  of  water 
vapor,  embodied  in  the  last  term  of  Eq.  (15).  By 
contrast,  any  initial  concave-up  curvature  (d2V/di2  > 
0)  in  cathode-supported  cells  must  entirely  be  due  to 
a nonohmic  activation  polarization.  The  analytical 
form  of  the  effect  of  concentration  polarization  in 
cathode-supported  cells  is  such  as  to  impart  a 
convex-up  curvature  (d2V/d/2  < 0)  over  the  entire 
range  of  current  densities  [6].  Thus,  if  the  activation 
polarization  is  ohmic,  the  curvature  of  V vs.  i in 
cathode-supported  cells  must  be  convex-up  (d2V/ 
d/'2  < 0)  over  the  entire  range  of  current  densities. 
Indeed,  this  has  been  observed  in  cathode-supported 
cells  at  1000°C.  In  anode-supported  cells,  however, 
the  initial  curvature  is  expected  to  be  concave-up 
(d2V7d( 2 > 0)  even  when  the  activation  polarization 
is  ohmic. 


4.  Summary 

The  present  work  has  shown  that  the  principal 
losses  in  SOFC  are  attributed  to  activation  and 
concentration  polarizations,  especially  when  the  elec- 
trolyte is  a thin  film  supported  on  an  electrode. 
Power  densities  as  high  as  1.9  fl/cm2  at  800°C  were 
obtained  with  a thin  film  (-10  p,m)  of  YSZ  de- 
posited on  a relatively  thick  (1  mm),  porous  Ni  + 
YSZ  anode.  The  activation  polarization  on  the 
cathode  side  can  be  substantially  decreased  by  using 
a composite  electrode  comprising  a mixture  of  an 
electrocatalyst,  such  as  LSM,  and  an  oxygen  ion 
conductor,  such  as  YSZ.  The  theoretical  analysis 
predicts  that  activation  polarization,  as  defined  in 
terms  of  an  effective  charge  transfer  resistance,  ReJtf, 
can  decrease  or  increase  with  an  increasing  cathode 
thickness  and  attain  a limiting  value.  For  the  set  of 
parameters  applicable  to  LSM  as  the  electrocatalyst 
and  YSZ  as  the  oxygen  ion  conductor,  the  ex- 
perimental results  show  that  Rf*  decreases  with  an 
increase  in  the  cathode  thickness.  Thus,  a composite 
electrode  of  LSM  + YSZ  is  very  effective  in  lowering 
the  activation  polarization.  It  is  known  that  activation 
polarization  at  the  Ni  + YSZ  anode  is  generally  much 
smaller  than  at  the  cathode.  This  suggests  that  at 
the  Ni  + YSZ  anode  must  be  very  small  and  that  a 
composite  electrode  of  Ni  + YSZ  is  very  effective  in 
lowering  the  activation  polarization  at  the  anode.  The 
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analysis  of  gas  transport  through  a relatively  thick 
anode  shows  that  the  V vs.  i traces  are  expected  to  be 
nonlinear,  in  accord  with  the  experimental  results. 
Further,  the  analysis  also  shows  that  the  initial 
concave-up  curvature  (dV/d(2>0)  in  anode-sup- 
ported cells  has  its  origin  in  both  activation  as  well 
as  concentration  polarizations.  The  present  results 
show  that  the  polarization  effects  in  SOFC  can  be 
minimized  through  a microstructural  optimization  of 
the  electrodes.  As  far  as  the  activation  polarization  is 
concerned,  the  use  of  composite  electrodes  spreads 
the  reaction  zone  from  the  electrolyte /electrode 
interface  some  distance  into  the  electrode.  In  this 
sense,  the  electrolyte/electrode  interface  is  a func- 
tionally diffuse  interface. 
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